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SUMMARY 


A transonic Wall-Interference Assessment/Correction (WIAC) 
procedure has been developed and verified for the 8- by 24-inch 
airfoil test section of the Langley 0.3-m Transonic Cryogenic 
Tunnel. This report is a user manual for the correction pro- 
cedure. It includes a listing of the computer procedure file as 
well as input for and results from a step-by-step sample case. 

INTRODUCTION 

A transonic computational procedure for assessing and 
correcting the wall interference for airfoils tested in the 8- 
by 24-inch slotted-wall test section of the Langley 0.3-m TCT 
has been developed. This procedure is based upon wall inter- 
ference assessment /correction (WIAC) ideas and codes of Kemp 
(refs. 1-5). His latest code, TWINTN4 (ref. 5), is the heart of 
the present procedure. This code includes an accounting for 
sidewall boundary-layer effects, derived from the work of 
Barnwell and Sewall (refs. 6-8), and a least-squares velocity 
matching criterion over the entire airfoil surface as proposed 
by Murman (ref. 9). 
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Adaptation of these ideas and codes into a semiautomated 
WIAC procedure (fig. 1) for the 0 . 3-m TCT was discussed in refer- 
ence 10. Validation of the WIAC procedure, using results from 
several tests in the 0 . 3-m TCT, was reported in reference 11. 
The reader should consult these latter two references for details 
about the WIAC procedure and results from it. Also included as 
part of this procedure is the GRUMFOIL (refs. 12-13) transonic 
full-potential equation airfoil analysis code with viscous 
interaction; this code provides an independent "free-air" check 
of the wall-interference corrections to Mach number and angle of 
attack which are determined by the TWINTN4 WIAC code. 



Figure 1. Schematic diagram of WIAC procedure for the 0.3m TCT. 
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Descriptions of the Langley 0.3-m TCT 8- by 24-inch airfoil 
test section and the advanced technology airfoil test program are 
included in references 14 to 19 as well as in the test data 
references, such as reference 20. The reader is also referred to 
these works for aspects related to the cryogenic tunnel concept, 
design, development, and test techniques. 

This user manual briefly describes the present WIAC proce- 
dure for the 8- by 24-inch airfoil test section of the 0.3-m 
TCT. A sample case is carried step-by-step through the WIAC 
procedure in Appendix B. 


DESCRIPTION OF PROCEDURE 

The WIAC procedure for data taken in the 0.3-m TCT airfoil 
test section is comprised of five FORTRAN programs : GETTEST, 
PRECOR, NOCORR, TWINTN4 , and GRUMFOIL . These programs begin 
with the retrieval of coefficient data from tapes made during the 
0.3-m TCT tests, then determine Mach number and angle-of -at tack 
corrections, and finally compare corrected data to an independent 
free-air calculation result. Figure 2 is a flow chart ‘which 
indicates the normal progression through this procedure. All of 
the routines are initiated interactively from a file containing a 
collection of procedures; this file is listed in Appendix A. 
To illustrate use of the numerous procedures, a sample case 
(test 169, run 13, point 8) will be followed through step-by-step 
from beginning to end in Appendix B. This sample case is for 
results from a test of the CAST 10-2/DOA 2 airfoil model. 
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Data to be corrected are retrieved from the coefficient data 
tapes made subsequent to the airfoil tests in the 0.3-m TCT. The 
code GETTEST accesses the tape, extracts the data required by 
the wall-interference assessment /correction code, and creates 
separate files for each run specified on the input file. This 
program is submitted by typing "-GETTEST ,, f llename , tapenum" where 
'filename' designates the input file and 'tapenum' designates the 
name of the coefficient data tape. When GETTEST has completed 
execution and created the requested files, then PRECOR can be 
executed as shown in figure 2. The code PRECOR is an interactive 
routine which uses DI-3000 graphics routines; therefore, it must 
be run on a graphics terminal. To initiate the program, the user 
should type "-PRECOR, , filename ,filename2 ,DD" where 'filename' is 
the name of one of the files created by GETTEST, 'filename2' is a 
file which describes the wake rake layout for the specific test 
and 'DD' is the device driver for the computer terminal being 
used. The file referred to as 'filename2' is a list of the wake 
rake tube numbers which were used in the drag survey. The naming 
convention for 'filename2' is PRCDttt where ttt is the three- 
digit test number. PRECOR writes the input files for each test 
point onto individual TWINTN4-type input files. The grid 
parameters for TWINTN4 are set for 20 grid points along the model 
chord located midway between the upper and lower limits of a grid 
depicting the 8- by 24-inch airfoil test section of the 0.3-m 
TCT. According to the file naming convention which has been 
adopted for this procedure, the files are called 'Attrrpp' where 
'A' designates a TWINTN4-type input file and 'ttrrpp' is the 
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concatenation of the last two digits of the test number, the two- 
digit run number, and the two-digit point number. For instance, 
for the sample case, the input file to TWINTN4 would be 
A691308. These files, while they are in the form to be read by 

TWINTN4, may also be used to run GRUMFOIL for a comparison of the 
uncorrected test data and the free-air results at the uncorrected 
test conditions. 

A program called NOCORR creates a GRUMFOIL input file from a 
TWINTN4-type input file. This step may be necessary as well as 
informative for several tests. The correction code TWINTN4 
requires data at the midpoint of each grid cell on the airfoil. 
For the grid developed for a 6-inch model, data must extend as 
far back as x/c = 0.975; however, several of the models tested 
in the 0.3-m TCT lacked a pressure orifice in that region. 
Therefore to run TWINTN4, the test data must be augmented. The 
free-air code GRUMFOIL, with a few additions, checks for data in 
the trailing-edge region; and, if necessary, computes a suggested 
value for the auxiliary data by extrapolating the measured 
pressures while assuming the same shape as the theoretical 
pressures at the uncorrected conditions. The suggested values 
can be added to the 'Attrrpp' file using the line editor. NOCORR 
is initiated by typing "-NOCORR ,, t t rrpp " which creates an input 
file called "Fttrrpp' and submits the GRUMFOIL job. 

The correction code TWINTN4 is submitted by typing 
"-TWINTN4 , , tt , rrpp . " The input files created by PRECOR direct 
TWINTN4 to run a unified four-wall correction and create an input 
file for GRUMFOIL called 'Bttrrpp.' The TWINTN4 input file can 
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be changed to run the sequential correction as described in 
reference 5. Also, TW1NTN4 creates a second GRUMFOIL input file 
called 'Ettrrpp' which reflects Only the Barnwell-Sewall type 
Mach number correction due to the sidewall boundary layer. The 
flow chart (fig. 2) indicates the second pass run of TWINTN4 with 
an updated upstream boundary condition. Finally, the free-air 

comparison is made by running GRUMFOIL which is submitted by 
typing "-GOWIMED , .Bttrrpp." 

DESCRIPTION OF INDIVIDUAL CODES 

GETTEST .- The code GETTEST calls a program called PF which 
was developed by the Systems Development Corporation to effi- 
ciently handle the data taken during 0.3-m TCT tests. The number 
of spanwise pressure orifices on the model and on the drag wake 
rake and the number of streamwise pressure orifices on the top 
and bottom tunnel walls and the upper and lower model surfaces 
must be specified on the input file. GETTEST creates a data file 
for each run of each of the tests which are requested from the 
coefficient data tape. 

PRECOR . - From a series of menu-driven inputs, PRECOR allows 
the user to preview pressures measured on the upper and lower 
model surfaces or on the top and bottom tunnel walls and, if 
necessary, to delete those at bad (leaking or blocked) pressure 
orifices . 

The first user input is namelist TAPE which consists of the 
following : 
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NMI(l), NMI ( 2 ) 


NWI(l), NWI ( 2 ) 


NCOR 


NSPA 


CHORD 

XTSHIFT 


DIMEN 


The next user 

NTEST 

NRUN 

NPOINT 

MPOINT 

IPLOT 


The number of streamwise pressure 
orifices on the upper and lower model 
surfaces, respectively: default 29,18; 
maximum 47,47. 

The number of streamwise pressure 
orifices on the top and bottom tunnel 
wall centerlines, respectively: 
default 26,28; maximum 47,47. 

The number of spanwise pressure orifices 
on the wake rake: default 5, maximum 6. 
The number of spanwise pressure orifices 
on the model: default 21, maximum 21. 

The model chord length: default 6.0. 

The chordwise location of the turntable 
pivot point with respect to the model 
leading edge: default 3.0. 

The scale factor to convert CHORD and 
XTSHIFT dimensions to inches; i.e., 1 for 
inches, 25.4 for millimeters: default 1. 

input is namelist IN which consists of : 

The test number. 

The run number. 

Array of point number(s). 

The number of test points selected: 
default 1, maximum 10. 

Plot trigger; 1 for plots, 0 to override 
plots: default 1. 
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If IPLOT is set to zero, no plots are generated; the editing 
capability is temporarily disabled and the TWINTN4 input files 
are created immediately. Otherwise, the first frame plotted 
shows the spanwise drag distribution and an oblique plot of the 
wind-tunnel model with pressure measurements off the model center 
line. On the following frame, plots of the model surface and 
tunnel wall pressure signatures are created and the user is given 
the opportunity (and instructions) to delete any number of points 
using the crosshair or similar graphics input device associated 
with the display terminal. The code then returns to the input 
namelist IN. A simple carriage return response exits the code, 
sends a table to a predetermined printer, sends the plot metafile 
to the Varian plotter, saves it as METrr , and writes separate 
input files for each of the selected test points. 

The naming convention used for the WIAC procedure data files 
is 'Attrrpp' where 'tt' denotes the last two digits of the test 
number, 'rr' denotes the two-digit run number, and 'pp' denotes 
the point number. This convention is carried throughout the WIAC 
procedure data files. The first letter identifies both the code 
which creates and that which uses the data file. The letter 'A' 
designates an input file for TWINTN4 created by PRECOR. The 
letter 'B' designates an input file for GRUMFOIL created by 
TWINTN4 which contains the Mach number and angle-of -attack 
corrections. The ' C' and 'D' letters designate plot files from 
TWINTN4 and GRUMFOIL, respectively. The letters 'E' and ' F' 
designate input files to GRUMFOIL with Mach corrections only and 
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no corrections, respectively. The 'E' file is created by TWINTN4 
and the ' F' file is created by NOCORR. 

TWINTN4 . - The code in the WIAC procedure which performs the 
actual wall-interference correction calculations is a version of 
TWINTN4 (ref. 5). Minor modifications and additions have been 
made to the input and output; these have no effect on the correc- 
tions. The standard input and output for the TWINTN4 WIAC code 
are described in the user manual, reference 5. The reader should 
consult reference 5 and earlier references 1-4 for information 
about this code. 

The input to TWINTN4 has been modified to include YM, the 
ordinate at each model pressure orifice. An input flag, IPRE, 
has been added which, when set to 1 during a sequential correc- 
tion (ISEQ*1), creates the input file for GRUMFOIL which takes 
into account the Barnwell-Sewall sidewall correction alone. 
Several additions have been made to the standard output. Two 
plots have been added : a comparison of the input model shape and 
the effective inviscid shape obtained as the inverse in-tunnel 
solution (measured pressures as boundary conditions) and a plot 
of the difference in camber distribution between the effective 
inviscid shape and the input model shape. These comparisons have 
also been added to the printed output. 

In the present WIAC procedure, TWINTN4 is run twice. The 
first run of TWINTN4 uses initial values for the upstream flow 
angles SLA and SLB which are calculated from the most upstream 
measured pressures. This first estimate neglects any angularity 
which might be present in the tunnel. A second estimate of 
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values for the upstream angularity inputs (SLA,SLB) is determined 


from the camber 

comparison of this 

first 

TWINTN4 run and 

is 

printed for 

use 

in the second run 

of TWINTN4. After 

the 

corrections 

have 

been completed, an 

input 

file is made 

f or 


GRUMFOIL . 

GRUMFOIL . - The code in this WIAC procedure which compares 
the corrected data with the independent free-air check is 
GRUMFOIL. The MCMJ-9 version (with modifications) is being used 
in the current procedure. The modified version adds another 
option to the IBL flag. Setting IBL = 2 indicates a WIAC run and 
the difference between the input model shape and the equivalent 
inviscid shape calculated in TWINTN4 is read from the standard 
input file. Since the standard small-disturbance style 
(Cartesian) grid in TWINTN4 is coarse and offers no clustering, 
the input which controls the mapping which occurs in GRUMFOIL is 
repeated. This allows a more accurate representation of the 
model shape to be read from a user-supplied file called ALDATtt 
where, as before, 'tt' are the last two digits of the test 
number. The repeated block of input begins with the title card 
for FSYM, NU , and NL (input line 15) as described in the GRUMFOIL 
input manual (Appendix C) and continues through the upper and 
lower surface coordinates. 

Upon completion of the computation, a comparison of the 
thickness and camber distributions is plotted. These are for the 
effective inviscid shape from TWINTN4 and the airfoil shape plus 
the boundary layer calculated by GRUMFOIL. This plot is an 
addition to the standard set of plots generated by GRUMFOIL which 
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includes a comparison of the free-air calculated model pressure 
coefficient distribution and the corrected tunnel data. 

SELECTION OF SAMPLE CASE 

The 0.3-m TCT test data sample case to be followed step by 
step through the WIAC procedure in Appendix B was selected in 

order to demonstrate not only the semiautomatic nature of the 
procedure, but also to show deviations which need special 
attention. For instance, the correction code requires pressure 
data between the x = 0.975 and x = 1.00 chord stations which were 
not available for the test case. An erroneous pressure 
measurement on the top wall needs to be discarded; this utilizes 
one of the abilities offered by the procedure. Due to severe 
interactions between the shock and the boundary layer, the 
equivalent inviscid shape becomes badly distorted; this often 
adversely affects the determination of the second-pass upstream 
boundary condition. This sample case is not a typical case in 
that it shows all of these problems. It should not be miscon- 
strued as a good example of the results obtained from the 

correction procedure; this case borders on being uncorrectable . 
This sample case is a point from the second test of the 6-inch- 
chord CAST 10-2/D0A 2 airfoil model (ref. 20) and is identified 
as test 169, run 13, point 8. The test data report is in 

preparation. This case is for the nominal tunnel test condi- 

tions M ~ 0.73, a ~ 4.0°, C L ~ 0.776. Additional results from 
application of the WIAC procedure to several tests of the 
CAST 10-2/DOA 2 are given in reference 11. 
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APPENDIX A 


Listing Of Procedure File For WIAC Method 

Following is a listing of the file which contains each of 
the procedures called by the user and those called by the other 
procedures. The user should change the delivery information on 
each of the JOB statements and supply the appropriate user 
number, password, and charge number to the VPS computer control 
statements and assign the appropriate printer designation. The 
procedures assume they reside on a file called WIPROC. The 
procedure file was written for the CDC NOS 1.4 operating system 
on the CYBER computers and the VSOS 2.1.5 operating system on the 
VPS-32 computer installed at the NASA Langley Research Center. 
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. PROC,PRECOR, J136022 ,PRCD136 ,DD4014. 

RETURN , EDFI L , S AVFI L , S AVPVF . 

GET, TAPE 1=J 136022. 

FILE , TAPE 1 , BT= I . 

GET , TAPE 7 “PRCD 136. 

BEGIN, GODI, WIPROC ,PRECOR,DD4014 . 

BEGIN , SND , WIPROC , TTABLE , L0=$B1244CG PREPROCESSOR TABLE $. 

REWIND, EDFIL,SAVFIL. 

XEDIT , TWININ , NH, I-EDFIL, L=0 . 

BEGIN, ,SAVFIL. 

REVERT. DATA PREPROCESSER J136022 

. PROC , GODI , FNAME , DD4 105 , NC=1 / 0 , L=0/LIST,MAP=N/ $ SB/ SB$ . 

IFE, (NC.EQ. 1) ,C0MP. 

RETURN, LGO. 

GET, FNAME. 

FTN5,I=FNAME,#L=L. 

ENDIF ,COMP. 

ATT AC 11 , DI3 000 , MFNODE , SSNODE , D IERFN , DD4 10 5 /UN=LIBRARY . 

LDSET ,LIB=DI3000 , #MAP=MAP. 

LOAD, MFNODE , SSNODE , DD4 105 , LGO . 

NOGO.DIGO. 

DIGO. 

REVERT. COMPILE ANDLOAD DI3000 

.PROC, SND, FNAME, EI=$ ,M=T, EI“101353$/$ ,M“R,EI=101333$ ,L0=B1244CG,EX=LP/IN. 

DELIVER.LO 

S END , FNAME_EI , DC=EX . 

REVERT. SEND _E I 

.PROC .TWINTN4 ,RN=1 92701 ,WITWMOD. 

BEGIN, SND, WIPROC, SENDEE, , ,EX. 

REVERT. TWINTN4: TEST, RUN, POINT RN .MODS 
.DATA, SENDEE. 

/JOB 

TWINTN4 , T2800 . B1244CG TWINTAN/GRUMFOIL 

/USER 

/CHARGE 

ENQUIRE, OP=B. 

GET , UTWINM2 . 

GET, WITWMOD. 

UPDATE , I “WITWMOD , P=UTWINM2 , F , C , L=0 . 

GE T , NMACFTN / UN =L IBRARY . 

ATTACH , LARCGOS/ UN=LIBRARY . 

FTN5 , I=COMPILE , B=CHANGE ,L=LTWIN . 

GET , TAPE5=A_RN . 

COPYSBF , TAPE 5 , OUTPUT . 

REWIND, TAPE5. 

LDSET , LIB=NMACFTN/LARCGOS , PRESETA=NGINDEF , MAP=SB/ TWMAP . 

LOAD, CHANGE. 
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NOGO , TEMP . 

TEMP (TAPE 5) 

REPLACE , TAPE3=B_RN. 

REPLACE, TAPE10=E__RN/NA. 

PLOT . VARI AN , FF 
REPLACE , SAVPLT=C_RN . 

EXIT. 

REPLACE, TWMAP/NA. 

REPLACE, LTWIN/NA. 

DAYF ILE , DAYT . 

REPLACE.DAYT. 

'/EOF 

. PROC,NOCORR,RN=3 51 005, ALDAT35, SET-0 5, DAT-135. 

IFE , F ILE ( BNOCORR , .NOT. AS) ,MOVON. 

GET, NOCORR. 

FTN5 , I-NOCORR, L-0 ,L0=0 , B-BNOCORR. 

ENDIF.MOVON. 

MAP (OFF) 

GET , TAPE5=A_RN . 

BNOCORR. 

REPLACE , TAPE3=F_RN . 

BEGIN , GOWIMED , WIPROC , F_RN, ALDAT3 5 , BWI ST , SET , DAT . 

REVERT. WIAC UNCORRECTED: TEST-RUN-PNT RN ALDAT #SET_SET 

. PROC , GOW IME D , B 1 9 0 1 0 1 , ALDAT 1 9 , B WI S T , SE T=1 , DAT-DAT . 

SUBMIT, SENDEE, B. 

REVERT. WIAC/ GRM9 J3190101 ALDAT BIN #SET_SET 
.DATA, SENDEE. 

/JOB 

STAR, T300, CM75000, STRHZ . B1244CG CLYDE GUMBERT 

/USER 

/CHARGE 

GET,B190101 . 

GET,TAPE3_SET=ALDAT19. 

TOVPS ( INPUT , C6UD-B1 90101 , TAPE3_SET , UN-#//// , PW- , AC-) 

EXIT. 

DAYF ILE, ERRDAY. 

REPLACE, ERRDAY. 

/EOR 

GRUMFL. 

USER , U =//// // , PA- , AC- . 

RES0URCE(TL=450 ,LP=4,WS=800, JCAT-MDBAT) 

DELIVER, B1244CG. 

ATTACH, BWIST, WAIT- Y, USER-////// . 

TV, 4. 

LOAD ( BWIST, CN-GO , #800 , CDF-//1000 ,LIST-MAPST, GRLPALL- ) 

COMMENT. TONOS( Z ,C6UD=MAPST , JCS="###N" , " PWORD" ,"CHNUM") 

ATTACH, B190101,TAPE3_SET. 

GO(TAPE5=B190101,TAPE7=PLW SET,TAPE6-0W_SET,TAPE3-TAPE3 SET,TAPE10-T10 SET) 
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TONOS(Z,C6UD=T10_SET, JCS-"###N" ,"PWORD" ."CHNUM") 
T0N0S(Z,C6UD=0W_SET, JCS-"###N" ,"PWORD" , "CHNUM") 
T0N0S(Z,C6UD=PLW_SET,JCS“"iW/N" ,"PWORD" , "CHNUM") 

PURGE, B190101,TAPE3_SET. 

TONOS ( INPUT) 

DAYF ILE , DW_SET . 

TONOS( Z ,C6UD=DW SET , JCS-"###N" , "PWORD" , "CHNUM") 

EXIT. 

PATTACH, UTILITY. 

PURGE ,TAPE3_SET , B1 9010 1 . 

DAYF ILE, DW_SET. 

TONOS ( Z , C6UD a *DW_SET , JC S=" ###N" , "PWORD" , "CHNUM”) 

IDUMP ( I =G0 , L=D Ul'lPST ) 

TONOS(Z,C6UD=T10_SET, JCS-"###N" , "PWORD" , "CHNUM") 

TONOS (Z,C6UD»DUMPST, JCS-"###N" ," PWORD" ," CHNUM") 

TONOS(Z , C6UD=OW_SET , JCS-"###N" , "PWORD" , " CHNUM") 
T0N0S(Z,C6UD-PLW_SET, JCS-"###N" ."PWORD" ."CHNUM") 

TONOS (INPUT) 

/EOR 

FILES, T200. B1244CG CLYDE GUMBERT 

/USER 

/CHARGE 

ATTACH , T1 0_SET/NA . 

GET, #DAT_DAT/NA. 

IF E , F I LE ( #DAT_DAT , AS ) .SKIP. 

SKIPEI , #DAT_DAT. 

ENDIF.SKIP. 

COPYE I , T1 0_SET , #DAT_DAT . 

PACK, #DAT_BAT. 

REPLACE, #DAT_DAT/NA. 

PURGE , T1 0_SET/ ST=LPF . 

GET.WIGRPP. 

ATTACH, PLW_SET. 

C OPY , PLW_S ET , TAPE 5 . 

REWIND, TAPE 5. 

FTN5,I=WIGRPP,L=*0. 

RETURN, WIGRPP. 

ATTACH(LARCGOS/UN=LIBRARY) 

MAP (OFF) 

LDSET(LIB=*LARCGOS , PRESETA-NGINF) 

LGO. 

REPLACE, SAVPLT=GRW_SET. 

ATTACH , OW_SET , DW_SET/NA . 

COPY , OW_SET .LIST. 

COPY, DW_SET, LIST. 

REWIND, LIST. 

BEGIN, SND.WI PROC .LIST. 

COMMENT. ROUTE, LI ST, DC =LP. 

PLOT. VARIAN,FF 

PURGE, OW SET.DW SET/ ST=LPF , NA. 
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EXIT. 

DAYF ILE , ERRDAY . . 

RE PLACE , ERRD AY . 

REWIND, OUTPUT. 

COPY, OUTPUT, OUT. 

BEGIN , SND , WIPROC , OUT , El . 

DAYFILE, ERRDAY. 

REPLACE, ERRDAY. 

REPLACE , SAVPLT=ERRPLT . 

/EOF 

. PROC , GE TTE ST , FNAME , TAPNUM . 

GET, FNAME. 

SUBMIT, SENDE, B. 

REVERT. GET TCT TAPE TAPNUM 

. DATA , SENDE 

/JOB 

/NOSEQ 

GETTEST , TIOOO , STRHC . B1244CG 

/USER 

/CHARGE 

LABEL( TAPE 1 , D=PE , LB=KL , F=I , NT , PO=RN, VSN-TAPNUM) 
REWIND ( TAPE 1) 

GET,PF/UN=474750C. 

PF(TAPE1,PL=10000) 

DAYF ILE , DAY 1 . 

REPLACE, DAY 1. 

EXIT. 

DAYFILE, ERRDAY. 

REPLACE, ERRDAY. 

/EOR 

/READ, FNAME 

/EOR 

/EOF 

. PROC , GRM9COM , UGRM9ST , BWIST , WIMDMOD , BE . 

SUBMIT, SENDE, B. 

REVERT . COMP-P=UGRM9 ST , B=BWIST , I-WIMDMOD ,OPT=BE 
.DATA, SENDE. 

/JOB 

STAR , T1 00 , STRHZ . B1244CG 

/USER 

/CHARGE 

GET , OLDPL«UGRM9ST . 

GET , WIMDMOD/NA. 

UPDATE , 1= WIMDMOD, F , C=BWIST_C . 

TOVPS ( INPUT , C 6UD= BWI S T_C , UN -### , PW= , AC- ) 

EXIT. 

DAYF ILE, ERRDAY. 

REPLACE, ERRDAY. 


TCT TAPE CREATION 


CLYDE GUMBERT 
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/EOR 

BWIST_C. 

USER( U-M# , PA= , AC-) 

RESOURCE (TL=*1 50 ,WS=660 ,LP**4 , JCAT“SMBAT ) 

DELIVER(B1244CG) 

TV, 8. 

ATTACH, BWIST_C. 

PURGE, BWIST. 

FORTRAN ( I=BWIST_C , B=BWIST/#400 ,L-COMPST , OPT-BE) 

TONOS(Z ,UUUD»BWIST, C6UD-C0MPST, JCS»"###N" , "PWORD" , "CUNUM" ) 

DEFINE, BWIST. 

FILES , PR=* , L=VPSLIST . 

SUMMARY. 

PURGE, BWIST_C. 

DAYF ILE , DCOMST . 

T0N0S(Z,C6UD=DC0MST,VPSLIST, JGS-"###N" , "PWORD" ,"CHNUM") 

EXIT. 

DAYFILE, DCOMST. 

T0N0S(Z,C6UD=DC0MST, JCS-"###N" , "PWORD” ,"CHNUM") 

/EOF 

. PROC , GLOT , DIMETA , CAL“0/ 1 , DI=0/ 1 ,E I“/EI . 

NOTE , SENDE , NR . \ / JOB\GOPLOT , T10 . B1244CG PLOTS 

NOTE, SENDE ,NR.\/USER\/CHARGE 
IFE ,DI.EQ. 1 ,DI3. 

IFE , CAL.EQ. 1 ,CALP. 

NOTE , SENDE , NR . \GET , DBIETA. \PLOT . #CAL , 1 1 // 0.3MM LEROY PEN // 

ELSE ,CALP. 

NOTE , SENDE , NR . \GET , DIMETA. \PLOT . VAR, F 
ENDIF,CALP. 

ELSE,DI3. 

IFE, CAL.EQ. 1,CALB. 

NOTE, SENDE , NR. \GET, SAVPLT=DIMETA.\PLOT. CALPOST, 11 // .3MM LEROY PEN // 

ELSE.CALB. 

NOTE , SENDE , NR. \GET, SAVPLT-DIMETA. \PLOT. VARIAN, FF 
ENDIF,CALB. 

ENDIF ,DI3 . 

NOTE , SENDE , NR . \EXIT . 

IFEjDI.EQ. 1,DI2. 

NOTE, SENDE, NR. \/EOR\S MF I DIMETA\D P FROM FIRST TO LAST\Q\/EOF 
ELSE,DI2. 

NOTE, SENDE, NR. \ /EOF 
ENDIF, DI2. 

PACK, SENDE. 

REWIND, SENDE. 

BEGIN, SND,WIPROC, SENDE, El, ,EX. 

REVERT. SEND PLOTS 
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APPENDIX B 


Sample Case 

This appendix shows the required user inputs such as input 
files for batch jobs or specific entries for interactive jobs for 
the sample case: test 169, run 13, point 8. The input file for 

procedure GETTEST was called GET169; it is listed in figure Bl. 
The coefficient data for test 169 were written on a tape named 
NJ0181. There were 26 and 18 chordwise pressure orifices on the 
model centerline upper and low surfaces, respectively. There 
were 6 probes on the drag rake and 21 spanwise pressure orifices 
on the model. There were 26 and 28 pressure orifices on the top 
and bottom wall, respectively. The 'FACILITY L' command 

specified the first letter of the created files to be an 'L.' 
The command to run GETTEST was "-GETTEST, WIPROC, GET169 ,NJ0181 . " 
Of the files created, the one which contains the data for the 
sample case is 'L169013.' This file can be accessed by the 

program PRECOR. 

The interactive terminal session during which PRECOR was 
executed is illustrated by prints of the screen, as shown in 
figures B2-B7. Figure B2 shows the command to initiate PRECOR, 
the printout of the default values in NAMELIST TAPE, and the 
updated values that were input. Figure B3 shows the next print 

which calls for input to NAMELIST IN and user response for this 
sample case. Since the variable, IPLOT, was let at the default 
value of 1, the code began to draw the plots. The first frame 

(fig. B4) shows the spanwise variation of the drag coefficient 

determined from the wake rake survey, and an oblique projection 
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of the wind-tunnel model above which pressure coefficients are 
plotted at all locations off the centerline. The experimental 
pressure distributions on the model and the walls were plotted as 
shown in figure B5. Instructions were then printed to allow 
editing of the data. 

There are no indications that the model pressures have any 
"bad" points; therefore, a space was entered for each of the 

first two prompts. However, the 14th orifice on the upper wall 
was classified a "leaker" and the value looks inconsistent with 
the rest; therefore, it was removed as shown in figure B6. 

PRECOR echoed the point which had been deleted. The first point 
on the upper wall is also bad; however, this point is expressly 
needed by the correction code so it must not simply be 

eliminated. Rather, after PRECOR has been run and the file 
A691308 has been created the value must be fixed using an 
editor. Experience gained from other tests in this tunnel 
indicates that, in general, the Cp variation over the first 

three orifices of the top wall is nearly linear. Therefore, the 


value 

to 

be used for 

the 

first point 

on 

the top 

wall can 

be 

linea 

r ly 

extrapolated 

from 

points two 

and 

three . 

Since th 

ere 

were 

no 

other points 1 

to be 

removed, a 

space 

bar was 

entered 

f or 


each of the two following prompts. The test, run, and point 
numbers were echoed as shown in figure B7. A carriage return is 
then required to continue to the next case. The prompt to input 
NAMELIST IN was answered by merely a carriage return. This 
indicated an end-of-file to the program, causing it to exit, to 
save the TWINTN4-type input file, A691308, and to send a data 
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table (Table Bl) and a plot file to their respective output 
devices. Information on Table Bl includes the identification of 
the test, run, and point; the test parameters; the spanwise drag 
and pressure distributions; and the two upstream wall-pressure 
measurements. Upon completion, the procedure echoes the message 
"REVERT. DATA PREPROCESSOR L169013" to the screen. 

The plots of the model pressure distribution shown in 
figures B5-B7 indicate that the data do not extend back to the 
x/c = 0.975 point designated by the cross on the x-axis. 

Therefore, before TWINTN4--the correction code — can be run, the 
data must be augmented. A means of extrapolating the experi- 
mental data has been added to GRUMFOIL, which can be run at this 
point by N0C0RR using the command "-N0C0RR, WIPROC , 69 1308 . " The 
extrapolated values of the pressure distribution are shown in 
figure B8 as they were printed by GRUMFOIL. These values were 
added to namelist variables YM, XM, and CPM and the values of 
NMI , the number of points, were increased to 27 and 19 on file 
A691308 using the line editor. 

To submit the correction code, enter the command 
"-TWINTN4, WIPROC, 69, 1308." Unless the input file, A691308, is 
changed, it will run a unified procedure assuming no upstream 
flow angularity. The comparison of the first-pass free-air model 
pressure distribution and the renormalized measurements is shown 
in figure B9. An additional result of the first run is the new 
suggested input value for SLA and SLB which can be changed on 
A691308. The suggested new values for SLA and SLB are derived 
from the comparison of the camberlines of the equivalent inviscid 
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shape and the model shape. One of the plots created by the 
modified version of TWINTN4 makes such a comparison. The differ- 
ence between the equivalent inviscid ordinates and the model 
ordinates is separated into a thickness part and a camber part. 
The camber part is shown for the sample case in figure BIO. For 
cases exhibiting strong shocks or regions of boundary-layer 
separation, the effective inviscid shape may become distorted. 
Figure Bll shows the camber comparison for point three of the 
same run for which the equivalent shape is not as distorted. For 
cases which exhibit distorted shapes such as the sample case, the 


suggested values 

of SLA and 

SLB 

may 

be wrong. 

However , the 

tunnel angularity 

appears to 

be 

(see 

ref. 10) 

nearly linearly 


dependent on the circulation and hence the lift. The second-pass 
values for SLA and SLB can be determined in such cases by 
extrapolating with respect to lift coefficient the values from 
lower lift cases of the same Mach number and Reynolds number as 
indicated in figure B12. The circles represent the slope of the 
camberline from the first pass; the diamond is the extrapolated 
value . 

The second run of TWINTN4 is submitted by again typing 
"-TWINTN4 ,WIPROC, 69 , 1308." The plots from the second run of 
TWINTN4 are shown in figures B13-B15. The first two figures are 
plots of the perturbation velocity at the tunnel centerline 
caused by all four tunnel walls (fig. B13) and the comparison of 
the free-air and renormalized model pressure distributions 
(fig. B14) as described in reference 5. The third figure, 
figure B15, is a direct comparison of the measured coordinates of 
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the model and the effective inviscid shape determined by 
the inverse solution for the tunnel flow using the wall and 
model pressures as boundary conditions. The last figure, 

figure B16, is a comparison of the camberline increment between 
the two shapes. 

The free-air check code, GRUMFOIL, is submitted to the 
VPS-32 by the command "-GOWIMED , WIPROC ,B69 1 308 , ALDAT69 . " The 

corrected Mach number and lift coefficient are specified and the 
drag and angle of attack are returned. The boundary-layer 
transition location is set close to the leading edge which seems 
to be appropriate for most cases at moderate Reynolds numbers. 
( See ref . 11.) 
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APPENDIX C 

Informal User Guide To GRUMFOIL Version MCMJ-9 

The methodology used in GRUMFOIL MCMJ-9 is described in 
reference 13. This appendix is a brief description of the input 
for the MCMJ-9 version of GRUMFOIL for which a formal users' 
guide is not available. 


Line 

1 

Title card (in 

A format) 

Line 

2 

header line 


Line 

3 

FNC, FNR, FMESH 

( F10 format ) 


FNC The number of cells around the airfoil for the 
crudest mesh. 

FNR The number of cells from the airfoil surface 
to farfield for the crudest mesh. 

FMESH The number of meshes to be run. 

Typical values are 40., 8., 3. 

Maximum dimensions allow 196 x 34 grid. 

Line 4 I CASE , I BL , I WCV , LMODE , IXP , IOPUT (11,14,515 format) 

ICASE Identifies run number. 

IBL Boundary-layer flag. 

0 for inviscid run. 

1 for viscous interaction. 

IWCV Wake curvature flag. 

0 includes wake curvature; 1 suppresses wake 
curvature . 

LMODE Kutta condition flag. 

-1 for inverse mode; c^ prescribed, returns a. 
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Line 5 


0 for direct mode; a prescribed, returns c^ . 

+1 for relaxed Kutta condition; both n and 
prescribed. Kutta condition not enforced. 

Run LM0DE=+1 only in inviscid mode (IBL=0); HMIN 
must also be specified. 

IXP Experimental data flag. 

0 no comparison with experimental data. 

1 experimental data included for comparison. 
IOPUT Output level flag. 

0 creates final plots only. 


7 draws plot for each grid, prints Mach 
chart and interacted boundary-layer 
solution after each update. 

Namelist &MXF 

(Note: VPS computer uses (&) instead of ($) for 

denoting and enclosing namelists.) 

WAL0,1,2 Under relaxat ion factor applied to a 

when LM0DE=-1 or to c^ when LMODE=0 on 
the (0,1,2) grid, 
defaults .5, .5, .5 

ITBLMO,l,2 Maximum number of update cycles on 
the (0,1,2) grid when IBL=l. 
defaults 22,22,30 
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NHDST 


NXNV 


LLO ,1,2 


LRO ,1,2 


NSMTH 


COVA 


COVL 


XPLT, YPLT 


Number of inviscid Iterations on the Oth 
grid before first boundary-layer update, 
default 2* TOT (from line 7) 

Number of iterations between a updates 
when LMODE=- 1 . 
default 4 

Number of points upstream of shock over 
which running average of surface velocity 
is taken for input to boundary layer. 
Negative value bypasses smoothing, 
defaults 1,1,1 

Number of points downstream of shock over 
which running average of surface velocity 
is taken for input to boundary layer. 
Negative value bypasses smoot ing. 
defaults 1 , 1 ,1 

Number of times smoothing is applied to 
surface velocity near shock waves, 
default 4 

Convergence criterion applied to angle of 
attack in inverse mode, 
default .005 

Convergence criterion applied to lift 
coefficient in direct mode, 
default . 005 

Dimensions of experimental-theoretical 
pressure comparison plot in inches. 
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defaults 5 . , 8 


Line 6 
Line 7 


Line 8 
Line 9 


CP1,CP2 Lower and upper limits of pressure 

comparison plot, 
default 1.2, -2.0 


Several other parameters are included in NAMELIST 
MXF; however, under all but the most strenuous 
conditions their default values will be adequate, 
header line 

TOT , COV , P 1 ,P2 ,P3 ,P4 ,P5 ,NMESH (F10 format) 

There will be a line 7 for each mesh (FMESH lines). 
Number of multigrid cycles. 

Convergence criterion applied to max 
residual . 

Relaxation factors; P4 is the primary 
relaxation factor typically ranging from 
.6 to .8; rarely change the others. 

Number of grids in multigrid cycle. 
Largest N for which 2 N is a factor of the 
smaller grid dimension. 

header line 

QC,VIS1 ,VIS2 ,VIS3 ( F10 format) 

QC Switch Mach number. 

VIS1 

VIS2 Artificial viscosity parameters. 

VIS3 


TOT 

COV 

P 1-5 
NMESH 
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Suggested values . 85 , 1 . 2 , 0 . , 2 . , respectively. 

Line 10 header line 

Line 11 AL1 , AL2 , STEP , FM1 ,FM2,DMACH (F10 format) 

AL1 Initial angle of attack in degrees. 

AL2 Final angle of attack. 

STEP Increment in angle of attack. 

FM1 Initial Mach number. 

FM2 Final Mach number. 

DMACH Mach number increment. 

In inverse mode AL1,AL2 and STEP may be left blank 
for single case run; if only AL1,FM1 are used, leave 
others blank. 

Line 12 header line 

Line 13 CLFT1 , CLFT2 ,DCLFT (F10 format) 

GLFT1 Initial lift coefficient. 

CLFT2 Final lift coefficient. 

DCLFT Lift-coefficient increment. 

For single case run, only CLFT1 is used, leave others 
blank. 

Line 15 FSYM , FNU , FNL (F10 format) 

FSYM |FSYM|>2 : input airfoil slopes. 

| FSYM | <2 {program calculates airfoil 
slopes . 

FSYM>0 {symmetric airfoil, input upper 
surface only. 

FSYMCO {asymmetric airfoil. 
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FNU 

Number of points to define upper 

surface . 



FNL 

Number of point to define lower 

surface . 

Line 

16 

header line 




Line 

17 

TRAIL. SLOPT 

(F10 format) 





TRAIL 

Trai ling -edge included 

angle in 

degrees . 




Program computes angle 

if TRAIL<0 and 




| FSYM | <2 . 





SLOPT 

Meanline slope at trailing edge, 

computed 




when |FSYM|<2. 



Line 

18 

header line 




Line 

19 

X.Y. SLOPE ( F10 format) FNU lines 





Upper surface ordinates and slopes 

from leading edge 



to trailing 

edge . 



Include 

lines 20 and 21 

only if FSYMCO. 



Line 

20 

header card 




Line 

21 

X.Y. SLOPE ( F 10 format) FNL lines 




Lower surface ordinates and slopes from leading edge 
to trailing edge. 

Line 22 Namelist &IINP 

Boundary-layer transition flag. 

1 Natural transition 

2 modes; KT=1 preferred 

3 by Melnik. 

4 Assigned transition mode; must 

input XTRANU ,L in &AINP. 

default 1 

Boundary-layer output flag. 


IOUT 


29 



-1 


Least output . 

0 

1 Most output, 

default -1 

Line 23 Namelist &A1NP 

REYNO Free-stream Reynolds number based on 

chord of 1.0. 

XTRANU , L Upper and lower surface assigned 

transition locations measured in x/c from 
the L.E. 
default . 10 

QEIT, B Upper and lower surface initial guess for 

transpiration velocity, order of 5 /c 
at trailing edge, 
defaults .004,, 002 

Increase to ,006 for REYN0<6 x 10^. 
Following lines (24-26) required if IXP^O (on line 4): 

Line 24 NXP (15 format) 

Number of experimental data points. 

Line 25 AMXP , ALXP , CLXP , CDXP (F10 format) 

AMXP Experimental Mach number. 

ALXP Experimental angle of attack. 

CLXP Experimental lift coefficient. 

CDXP Experimental drag coefficient . 

Line 26 XXP ,CPXP (F10 format ) 

XXP x-coordinates of experimental data. 

CPXP Experimental pressure coefficient. 
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To stack multiple case, start again at line 1. 

To end : 

Line 27 Title (END OF CALCULATION) 

Line 28 header line 

Line 29 0. (FNC=0. stops execution and prints closing 

banner.) 

The files needed to run GRUMFOIL MCMJ-9 are accessed on 
UN=350607N. For the CYBER VPS-32 version of GRUMFOIL MCMJ-9 
UGRM9ST Update OLDPL file 

BWIST VPS-32 object file 

WIGRPP Plotting routine 
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DATA TAPE FROM .3 M TCT TEST169 , RUN 13 , POINT 8 
TUNNEL CONDITIONS : REF MACH NO = .730 , REYNOLDS NO =10.0MILLI0N , ALPHA ■ 4.016 
CN = .7776 , CC =-.0083 , CMQC =-.0646 , CL = .7763 , CD = .0239 , CDP = .0461 

DRAG RAKE SURVEY : 

.0239 .0186 .0131 .0099 .0193 .0840 


SPANWISE PRESSURE SURVEY : 


X 

.505 

.505 

.505 

.505 

.505 

.505 

.505 

Y 

2.953 

2.559 

2.165 

1.772 

1.378 

.984 

.591 

Z 

.032 

.032 

.032 

.032 

.032 

.032 

.032 

CP 

-.619 

-.640 

-.635 

-.617 

-.652 

-.647 

-.646 


X 

.505 

.505 

.505 

.505 

.505 

.505 

.505 

Y 

.197 

-.197 

-.591 

-.984 

-1.378 

-1.772 

-2.165 

Z 

.032 

.032 

.032 

.032 

.032 

.032 

.032 

CP 

-.634 

-.613 

-.589 

-.569 

-.566 

-.581 

-.579 


X 

.505 

.505 

.064 

.145 

.225 

.305 

.385 

Y 

-2.559 

-2.953 

1.795 

1.795 

1.795 

1.795 

1.795 

Z 

.032 

.032 

.037 

.046 

.048 

.046 

.042 

CP 

-.056 

-.568 

-1.073 

-1.328 

-1.419 

-1.388 

-1.328 


TUNNEL WALL 
XT -26.500 

PRESSURE 

-24.500 

-22.500 

-20.500 

-18.500 

X -3.917 

-3.583 

-3.250 

-2.917 

-2.583 

CPU -.1268 

-.0111 

.0019 

.0033 

.0032 

CPL -.0350 

-.0079 

-.0013 

.0021 

.0115 


TABLE BI. - AUXILIARY DATA TABLE 
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FACILITY L 

TEST RUN POINT ALPHD M,INF RINF CN CC CMQC 
CPUSOl THRU CPUS26 USXOl THRU USX26 
USYOl THRU USY26 USZOl THRU USZ26 
CPLSOl THRU CPLS18 LSXOl THRU LSX18 
LSYOl THRU LSY18 LSZOl THRU LSZ18 

CPSPANOl THRU CPSPAN21 SPXOl THRU SPX21 SPYOl THRU SPY21 
SPZOl THRU SPZ21 CDCOR1 THRU CDCOR6 SBXOl THRU SBX28 
CPSBOI THRU CPSB28 STXOl THRU STX26 CPSTOl THRU CPST26 


Figure B 1 . Listing of GET169'. 


-PRECOR, , 

L169013 

-$TAPE 

ONMI 

s 

£9, 18 

0NU I 

S 

£6, £8 

0NCOR 

= 

5, 

0NSPA 

“ 

SI, 

0CHORD 

ss 

. 6E+01 

0XTSHIFT 

s 

. 3E+0 1 

0DIMEN 

0$END 

= 

. 1E+01 


’ *TAPE NMI=£6,NC0R=6, * 


Figure B2. Interactive execution of PRECOR 
- model/tunnel specifications. 


INPUT TEST NO , RUN NO , NO POINTS AND POINT NOS 
IN NAMELIST IN ... I.E. 

*IN NRUN= , NTEST = ,MPOINT= ,NPOINT= , , $ 

"> *IN NTEST=169,NRUN=13,NP0INT=8, * 


Figure B3. Interactive execution of PRECOR 
- choice of test, run and point number. 


# 
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Figure B4 . Interactive execution of PRECOR - spanwise drag distribution 
and upper surface pressure measurements. 




LOCATE CURSOR AND HIT 2ER0 TO REMOVE POINT 
FROM UPPER SURFACE(SQUARES) 

Z FOR NEXT CASE OR 

ANY OTHER KEY TO CONTINUE 


w 






Figure B5 . Interactive execution of PRECOR - 
pressure signatures. 


model and tunnel wall 



LOCATE CURSOR AND HIT 2ER0 TO REMOVE POINT 
FROM UPPER SURFACE (SQUARES) 



Figure B6. Interactive execution of PRECOR 


edit "bad" wall pressure 


LOCATE CURSOR AND HIT ZERO TO REMOVE POINT 
FROM UPPER SURFACE (SQUARES) 

Z FOR NEXT CASE OR 

ANY OTHER KEY TO CONTINUE 



Figure B7 


Interactive execution of PRECOR 


end of case 



.9614 


.946078 

•003b3y 

.bb54Bi 


.9705 


.954849 

.001508 

.658005 

.911652 

.9788 


.962929 

-.000457 

.651053 

.902992 

.9864 


.970290 

-.002249 

.644497 

.894801 

.9932 


.976905 

-.003861 

.638428 

.887212 

.9992 


.982750 

-.005285 

.631750 

.878681 

1.0044 


.987798 

-.006516 

.625411 

.870790 

1.0087 


.992022 

-.007546 

.617310 

.860096 

1.0122 


.995395 

-.008369 

.610022 

.851392 

1.0148 


.997882 

-.008976 

.600323 

.837931 

1.0164 


.999447 

-.009357 

.590040 

.826050 

1.0169 


1.000000 

-.009492 

.576647 

.804632 

X- .57836 

AT 

THE SHOCK 




CP- -.88559 

UPSTREAM 




CP- -.45953 

DOWNSTREAM, RANKINE-HUGONIOT 



CP- -.44259 

DOWNSTREAM, ISENTROPIC 



NXP.NC02: 44 

80 




XXP.CPXPI : 

4. 

7388 4, 

.7454 .1718 

.0626 


X,CP : 1.0000 

.9994 

.2110 .3313 


INPUT POINTS : 

* 

7 12 16 

NEIGHBOR X : 

.9478 .8977 



AUXILIARY POINT ADDED TO LOWER SURFACE 
EXPERIMENTAL DATA / . 

X , Y , CP .9795 -.0486 .1792 

INPUT POINTS : 155 145 141 NEIGHBOR X : 1.0000 .9001 

AUXILIARY POINT ADDED TO UPPER SURFACE 
EXPERIMENTAL DATA 

X,Y ,CP : .9812 -.0391 .2527 <p 


Figure B8. ^ Extrapolation of Cp near the trailing edge - 
GRUMFOIL output, uncorrected run. 


* 
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Figure 


4 

Figure BIO. 



0 .8 .4 .6 .8 1.0 


X/C 

THINTAN INPUT FROM THE NASA LANGLEY .3 METER CRYO TUNNEL 
MINF= .7113H CL= .73917 SIGCP= .10317 
A42AGLR 85/08/22. 07.66.52. 

B9 . TWINTN4 output, first pass C p comparison. 



.3 METER CRYO TUNNEL TEST 163 RUN 13 POINT 8 

Camber comparison after first pass - sample case. 
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THINTAN INPUT FROM THE NASH LANGLEY .3 METER CRYO TUNNEL 
SIGU= .00776 SIGVs .00331 SIGR= .00663 

mencon 66/02/14. 07.27.44. 

Figure B13. TWINTN4 output, second pass perturbation velocities 



THINTAN INPUT FROM THE NASA LANGLEY .3 METER CRYO TUNNEL 
MINF= .72546 CL= .78277 8ICCT= .04880 
A42ACQM 86/02/14. 07.27.44. 

Figure B14. i TWINTN4 output, second pass C D comparison. 




x/e 

ALPHA = 4.016 (DEGREES) ALP CORR = -.027 (RADIANS) 


Figure B15 


TWINTN4 output, second pass shape comparison. 
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.3 METER CRYO TUNNEL TEST 163 RUN 13 POINT 8 

Figure B16. ^ TWINTN4 output, second pass camber comparison. 
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